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Abstract 

A recent composite-dark-matter scenario assumes that the domi- 
nant fraction of dark matter consists of O-helium (OHe) dark atoms, 
in which a lepton-like doubly charged particle O is bound with a 
primordial helium nucleus. It liberates the physics of dark matter 
from unknown features of new physics, but it demands a deep under- 
standing of the details of known nuclear and atomic physics, which 
are still unclear. Here, we consider in detail the physics of the binding 
of OHe to various nuclei of interest for direct dark matter searches. 
We show that standard quantum mechanics leads to bound states in 
the keV region, but does not seem to provide a simple mechanism 
that stabilizes them. The crucial role of a barrier in the OHe-nucleus 
potential is confirmed for such a stabilization. 

1 Introduction 

Direct searches for dark matter have produced surprising results. Since the 
DAMA collaboration observed a signal, several other collaborations seem 
to confirm an observation, while others clearly rule out any detection. We 
summarize the situation in Table Q] and the current experimental situation 
is reviewed in [T]. This apparent contradiction comes from the analysis of 
the data under the assumption that nuclear recoil is the source of the signal. 
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Starting from 2006 it was proposed [21 [3] that the signal may be due 
to a different source: if dark matter has weakly bound states with normal 
matter, the observations could come from radiative capture of thermalized 
dark matter, and could depend on the detector composition and temper- 
ature. This scenario comes naturally from the consideration of composite 
dark matter. Indeed, one can imagine that dark matter is the result of 
the existence of heavy negatively charged particles that bind to primordial 
nuclei. 

Cosmological considerations imply that such candidates for dark mat- 
ter should consist of negatively doubly-charged heavy (~ 1 TeV) parti- 
cles, which we call O , coupled to primordial helium. Lepton-like techni- 
baryons, technileptons, AC-leptons or clusters of three heavy anti-U-quarks 
of 4th or 5th generation with strongly suppressed hadronic interactions are 
examples of such O particles (see [21 [3] for a review and for references). 

The cosmological and astrophysical effects of such composite dark mat- 
ter (dark atoms of OHe) are dominantly related to the helium shell of OHe 
and involve only one parameter of new physics — the mass of O . The 
positive results of the DAMA/Nal and DAM A/LIBRA experiments are ex- 
plained by annual modulations of the rate of radiative capture of OHe by 
sodium nuclei. Such radiative capture is possible only for intermediate-mass 
nuclei: this explains the negative results of the XENON100 experiment. 
The rate of this capture is proportional to the temperature: this leads to 
a suppression of this effect in cryogenic detectors, such as CDMS. OHe 
collisions in the central part of the Galaxy lead to OHe excitations, and 
de-excitations with pair production in E0 transitions can explain the excess 
of the positron-annihilation line, observed by INTEGRAL in the galactic 
bulge. 

These astroparticle data can be fitted, avoiding many astrophysical un- 
certainties of WIMP models, for a mass of O ~ 1 TeV, which stimulates 
searches for stable doubly charged lepton-like particles at the LHC as a test 
of the composite-dark-matter scenario. The problem with OHe dark matter 
is that its constituents may interact too much with normal matter. OHe 
is neutral, but a priori it has an unshielded nuclear attraction to matter 
nuclei. To avoid the problem, it was assumed that the effective potential 
between OHe and a normal nucleus would have a barrier, preventing He 
and/or O from falling into the nucleus, allowing only one bound state, 
and diminishing considerably the interactions of OHe. Under these condi- 
tions elastic collisions dominate in OHe interactions with matter, which is 
important for many aspects of the OHe scenario. 

In this paper, we show that indeed such a barrier is needed to make the 
model work, and we try to establish its existence through several methods. 
In the first section, we review the classical description of the problem [3] 
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Xe 


124-134 


54 


cryogenic 





Table 1: Results of various dark matter searches and composition of the 
detectors. 



and show that in fact it does not lead to a repulsive force. In section 2, 
we explore the spectrum of the bound states of OHe. We show that, if one 
considers only the screened Coulomb force, then bound states exist only for 
light nuclei, whereas if we consider a polarization of OHe due to a second- 
order Stark effect, then most nuclei have keV bound states. In the last 
section, we check that the description of the Stark effect that we used is 
reasonable via a perturbative calculation at large distances, but it is not 
reliable when the nucleus comes close to OHe, as one would then need to 
take into account a strong and inhomogeneous deformation of the ground 
state by the common effect of Coulomb and nuclear force. 



2 Classical model 

To study the polarization of the OHe atom under the influence of an ap- 
proaching A nucleus, we can first treat OHe as a classical structure and 
neglect the effects of O and nucleus motion. The polarization of OHe is 
then fixed by the equilibrium of forces acting on the He nucleus. For every 
position of the A nucleus, we can work on the O-A axis, in the rest frame 
of the O particle, as shown in Fig. [TJ 

We take the He and A nuclei as uniformly charged spheres of radii Rue 
and Ra and of charges Zn e = 2 and Za- We also assume that O is 
point-like. We then obtain the electrostatic potential for the interactions 
with O : 

Voa(za) = - ^ Q " , for za > Ra 

ZA 

-Z A Z a ( n z\ 



2R A V WJ^° IZA<Ra (1) 
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Figure 1: One-dimensional OHe atom polarized along the O —A axis, 
denoted z. zq is the distance between O and He and za is the distance of 
the A nucleus along the 2;— axis. 

for the interaction between the O and the A nucleus, and a similar ex- 
pression for VoHe(zHe)- The potential between the two nuclei has both 
electrostatic and nuclear contributions. In the former, we neglect the He 
size, and for the latter we use an experimental parametrisation of the a- 
nucleus potential from scattering experiments [9]: 

t r , x Z He Z A a -V Q 

V H eA{ZA ~ Z ) = 1 + -1 , nTTZ ZITZr w a . for \ Z A ~ ^0 > Ra 

\za — zq\ 1 + e^ ZA 20 H *)/ a 
Z He Z A a ( {z a -zq) 2 \ -V 



+ 



2R A \ R\ l + e (\ z A-zo\-R*)/a> 

for \z A - zo\ < Ra- (2) 

The nuclear interaction is represented in a Woods-Saxon form, with param- 
eters V = 30 MeV, a = 0.5 fm and R*{fm) = 1.35 x A 1 / 3 + 1.3 (fm). 

The equilibrium position zq of the He nucleus will be at the minimum 
of the potential VmA + VoHe an d will depend on z A . At that point, the 
Coulomb force balances the nuclear force: 

f He + F^HeA = ■ (3) 

When the equilibrium position is determined, the OHe- A potential is ob- 
tained by adding the dipole potential to the Woods-Saxon one 

VoHeA{zA) = V dip (z A ) + V WS {Z A - Zq), (4) 
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Figure 2: Polarization zq of the one-dimensional OHe atom as a function of 
the distance za of an approaching sodium nucleus. The break in the curve 
corresponds to the fact that for a distance of about 9 fm, He falls into the 
nuclear potential of A. 



is the dipole potential of the polarized OHe atom. 

Fig. [2] shows the polarization zq of the OHe atom as a function of the 
position za along the z— axis for an approaching sodium nucleus with Za = 
11. We see that it is negative at large distance, giving rise to an attractive 
dipole potential and that the nuclear force starts to reverse the dipole when 
the nucleus gets closer to O . 

This situation corresponds to a repulsive dipole that acts against the 
nuclear force, but it can be seen in Fig. [3] that this repulsive force is not 
sufficient to overcome the nuclear force between the two nuclei and to give 
rise to a repulsive global potential. 

When za ^$ 7.5 fm, Equation ([3J projected along the z— axis loses its 
initial solution and another one remains, that is located in the nuclear well, 
giving rise to a jump in the polarization and therefore in the total potential. 
Similar results and pictures can be obtained for other nuclei. 

Hence we see that, classically, no repulsive potential appears, even if the 
electrostatic force of OHe repels the A nucleus. In fact one can argue that 
this is a generic classical result which does not depend on the details of the 
calculation. If the configuration of the 3 objects is He-O-A, then clearly the 
force is attractive. If the configuration is O-He-A, then that means that the 
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Figure 3: Total OHe-A potential for sodium. 

nuclear force on He is larger than the electrostatic force from A. Again, a 
net attraction between OHe and A results. 

To settle this classical picture of permanent attraction in the OHe- 
nucleus system, a proper quantum treatment of the problem is needed. On 
the one hand, we shall see from the following discussion that simple semiclas- 
sical and perturbative descriptions cannot solve the problem of permanent 
attraction in the OHe-nucleus system. On the other hand, a crucial point 
may be missing in such treatments: the correct description of the nuclear 
effects when the nucleus is close to helium and when neither semiclassical 
nor perturbative approaches are valid. 

3 Semiclassical model 

The quantum problem involves very different scales: the OHe binding is of 
the order of one MeV, and we are looking for bound states of about one 
keV. To obtain both in the same framework would imply a solution of the 
3-body problem at better than one per thousand, which is clearly very hard. 

Fortunately, for the very excited bound states, one can use a simplified 
method. For these states, the OHe atom will not dissociate, so we can treat 
that system as a whole, allowing a small polarisation in the A direction. 
Furthermore, the interaction potential OHe-nucleus can be taken as radial, 
as the polarisaion of OHe will be in the A direction. Hence we can use 
spherical coordinates, with the O fixed at the origin and the position 
of the center of the nucleus A. We know in this case that the solutions of the 
Schrodinger equation take the form ^(^(r,^) = '^Myj- Y l m (9, <p) where 
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YJ m (#,<^) are the spherical harmonics and where the radial part Uk,i(r) has 
to satisfy the radial Schrodinger equation 



d 2 



Uk,i(r) 



+ 2m A E kjl -V(r) 




Uk,i(r) = 



(6) 



dr 2 



where I is the relative angular momentum, E^i is the total energy in the 
center-of-mass (O ) frame, and V(r) is the sum of the nuclear and of the 
electrostatic potentials between OHe and A. 

The next simplification comes from the fact that one is looking for 
weakly bound states, for which the WKB method applies, and consider- 
ably simplifies the solution. Finally, we further simplify the problem by 
approximating the He wave function in the OHe bound state by a Is hy- 
drogenoid wave function. 

The A nucleus is seen as a uniformly charged sphere of charge Za and 
of radius ii J 4(fm)= l.SSA 1 / 3 ^, where A is the number of nucleons in the 
nucleus. Its mass rriA is corrected by the nuclear binding energy B given 
by the Bethe-Weizsacker formula: m A = ^m,, + N A Tn n — B, where Na is 
the number of neutrons, and m p and m n are the masses of the proton and 
of the neutron respectively. 

The interactions between the OHe atom and the nucleus take two forms: 
nuclear attraction between the helium and the A nucleus at distances r < 
Ra and electrostatic interaction due to the electrical charges of the compo- 
nents at distances r > Ra- 

Out of the nuclear region, the electrostatic interaction is dominant and 
can be separated into two contributions : 1) the electrostatic interaction 
between the spherical charge distribution of the OHe atom in its ground 
state and the spherical charge distribution of the nucleus; 2) the electrostatic 
interaction between the polarized OHe atom and the nucleus due to the 
Stark effect. Therefore, we can write 



where Vcoul corresponds to Coulomb attraction between O screened by 
the helium charge distribution and A, and V stark represents the interaction 
term of the charged nucleus and the dipole. 

Outside the nucleus, i.e. for r ^ Ra, we find for the Coulomb term 



V E i 



Vcoul + Vstark 



(7) 



Vcoui(r) 



3 / -Z Z A a \ _2r/ 




(8) 
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with the Bohr radius of the OHe atom vq = 1/ 1 {mn e ZoZn e OL) ~ 1.81 fm and 
p = Ra/tq. This expression can be considered as an improvement of the 
form from |10] where the nucleus was assumed to be a point-like particle. 

For the Stark potential, we use the formula for the quadratic effect in a 
constant electric field |llj . taken to be the field of the nucleus at the position 
of O. The dipole moment of the OHe atom in its perturbed ground state 
can then be written: 

< ql >= (9) 

so that, for r Ra, 

Vstark = -qt.t = -\rlE* = -\rl^. (10) 

Expressions (|5J) and (|10p are valid when the nuclear effects are negligible. 

In the nuclear region, we take a trapezoidal nuclear well, which will 
simplify the WKB solution: 

V nuc i = -Vq for r < R* 

= v Eiec(R*+2a)+v (r _ R ^ _ y Q for ^ < r < ^ + 2a (11) 

= for r > i?* + 2a 

characterized by its depth Vq and its diffuseness parameter a representing 
the region of r in which it goes linearly from — Vq to VEiec{R* + o). From 
diffusion experiments of a particles on nuclei [9], one gets Vo ~ 30 MeV for 
nuclei with Za < 25 and Vo = 45 MeV for Za > 25, as well as a = 0.5 fm. 
In the following, we shall avoid the transition region Za G [21,29], which 
does not contain any nucleus used for direct dark matter detection. The 
nuclear radius parameter R* = {Ra + 1-3 + tq) fm, where Ra + 1-3 fm is 
taken from [9] to take the finite size of the alpha particles into account. 
Fig. [J] shows the form of the potential between OHe and Na. 

To find the spectrum corresponding to the potential of Fig. we use 
the approximate WKB solutions, which, once applied to each region, give 
a quantization condition for the energy. 

For I = 0, we obtain 




where p = J k(r)dr, b is the turning point such that E = V(b) and k(r) = 
^2m A (E - V(r)). 

At I 7^ 0, we know from |12j that the behaviour of the effective potential 
V e fj(r) = V(r) + 2m + r^ a ^ ^ e or isi n requires to modify the WKB method 
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Figure 4: Shape of the interaction potential OHe-nucleus, according to 
Eqs. [7] and [LTJ 



in the radial equation : V e ff(r) — > V e ff(r) = V(r) + \ m 2 Ar -i ■ Therefore, the 
quantization condition becomes 



where p = J a k(r)dr, a = Jq n(r)dr, a and b are the turning points such 



that E = V eff {a) = V eff (b) and k(r) = J2m A ( E - V eff (r) ) , k{t) = 



3.1 Spectra from a screened Coulomb potential 

In [3], the spectrum was considered for a screened Coulomb potential at 
long distance, as in Eq. [SI We reanalyse this question with our WKB 
formalism. For small nuclei Za < 20, we first fix the exact value of Vq to 
obtain the highest level at -3 keV for 23 Na from DAM A for / = 0. We 
obtain Vo = 31.9 MeV, in good agreement with [9] and use this value for all 
nuclei with Za < 20. The spectrum of the OHe- 23 Na system is shown in 



We see a rich spectrum with many levels in the MeV region, corre- 
sponding to nuclear levels, for which the WKB approximation may not be 
appropriate. The only level in the keV region is at I = 0. It can be con- 
sidered as being due to the presence of the electrostatic potential. It is 




2e 2<J cos p + sin p = 



(13) 
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Figure 5: Spectrum of the OHe- 23 Na system at different values of I. The 
energies are in absolute value. Vq = 31.9 MeV and a = 0.5 fm. 



remarkable that other nuclei such that Za < 20 do not have keV bound 
states. Fig. [6] shows the highest level at I = for the most stable nuclei for 
Za going from 1 to 20. It turns out that only 23 Na at Z^ a = 11 has a level 
in the keV region. 

For large nuclei Za > 30, the data indicate that the nuclear well is 
deeper. In this case, we take as a reference germanium Zee = 32, A = 74 
from CoGeNT, for which we find a highest level at I = in the keV region 
for Vq = 45 MeV, which is precisely the central value from [9] for larger 
nuclei. This value is used for all nuclei with Za > 30. Fig. [7] represents the 
spectrum of the OHe- 74 Ge system. It is of the same kind as for 23 Na, with 
only one level in the keV region. 

The second column of Table [2] gives the highest-energy level at I = for 
the large stable nuclei involved in the experiments of interest. According 
to this model, iodine and thallium from DAMA each admit one level in the 
keV region, while xenon from XENON100 doesn't. 
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Figure 6: Highest-energy level at I = for stable nuclei from Za = 1 to 
Za = 20. The energies are in absolute value. Vq = 31.9 MeV and a = 0.5 fm. 
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Figure 7: Spectrum of the OHe- Ge system at different values of I. The 
energies are in absolute value. Vq = 45 MeV and a = 0.5 fm 
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Nuclei 


E(keV) for a screened 
Coulomb potential 


E(keV) for a screened Coulomb 
potential added to a Stark potential 


r4 Ge 


3.88 


1.16 


127j 


0.500 


2.31 


132 Xe 


540. 


2.33 


184 W 


350. 


1.86 


201rp| 


15.6 


52.7 



Table 2: Highest-energy level at I = for some heavy stable nuclei from the 
experiments of interest, when Veu c = Vcoul (second column), or Veuc = 
Vcoul + Vstark (third column). The energies are in absolute value. Vq = 
45 MeV and a = 0.5 fm. 



3.1.1 Spectra from a screened Coulomb potential and a Stark 
potential 

The results can de discussed in the same way when Vstark + Vcoul is used 
in the calculations, and the values of Vo are identical to the central exper- 
imental values, i.e. 30 and 45 MeV, for small and large nuclei respectively. 
Fig. [8] illustrates the results in the particular case of the OHe- 23 Na system. 

The major difference lies in the fact that, in this case, the levels in the 
keV region are obtained more easily, with sometimes several keV levels for 
the same nucleus, especially for large nuclei. The reason lies in the shape 
of Veicci that is deeper and less steep when Vstark is used. Fig. [9J as well as 
third column of Table [2] show that most nuclei now have keV bound states. 
Hence the inclusion of the Stark potential seems to destroy the previous 
interpretation of the data, which relied on Na and Ge to be very special 
nuclei. 



4 Perturbative analysis 

The three-body OHe-nucleus bound-state problem can be simplified in an- 
other way, by noting that helium is much lighter than the A nuclei. Given 
this, one can simplify the total hamiltonian of the system, written in the 
reference frame of the O particle, and choosing the z axis in the direction 
of A to: 

H « -— !— Ai + VoHe (n) + V OA (R) + V H eA (h 2 ) (14) 
2niHe 

in which the kinetic energy term of the A nucleus has been neglected and 
where r \ is the position of the He nucleus, R is the distance of the A nucleus 
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Figure 8: Spectrum of the OHe- 23 Na system at different values of I, when 
^Elec = Vstark + Vcoul- The energies are in absolute value. Vo = 30 MeV 
and a = 0.5 fm. 
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Figure 9: Highest level at I = for stable nuclei from Za = 1 to = 20 
due to an electromagnetic potential Vstark + Vcoul- The energies are in 
absolute value. Vq = 30 MeV and a = 0.5 fm. 
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on the positive part of the z— axis and r\% is the distance between He and 
the A nucleus. Vjj stand for the interaction potential between I and J. We 
are thus left with the one-body problem of He in a total potential depending 
on the parameter R. We shall consider here the contribution of the external 
A nucleus as a perturbation to the OHe atom. The hamiltonian (|14|) can 
be rewritten as the sum of an unperturbed part Hq and a perturbation W: 

H = H + W (15) 

where Hq = — 2 m H + VoHe i r i) corresponds to the isolated OHe atom 
and where W = VoA (R) + Vile A (^12) is due to the presence of the external 
nucleus. We use here the following two-body interaction potentials : 

VoHe{ri) = 

Vqa(R) = ,R>R A 



R 

ZpZACt [ q R 2 
2R A 



3 - #r , R < Ra 



V HeA (r 12 ) = + l+e{ ~ 2 \ A)/a ,r 12 >R A 

= ^ ( 3 - Rt) + >^<^ 

where, in order to simplify the calculation, the OHe atom is treated as a 
hydrogenoid system, and where the A nucleus is seen as a sphere of radius 
itU(fm)= 1.35 x A 1 / 3 [9] and charge Z A in the definitions of potentials Voa 
and VueA that are therefore point (O or He ++ ) - sphere (A) interaction 
potentials. A Woods-Saxon potential of parameters Vq and a has been 
added in VueA to take the nuclear interaction of both nuclei into account. 

We are studying the perturbed ground state Eq(R) of the OHe atom un- 
der the influence of the external perturbation W(R). The perturbed energy 
is therefore an approximation of the energy of the total O-He-A system, 
described by the hamiltonian (|15p . If this energy presents a minimum for 
some Rh, then the system will tend to this configuration to minimize its 
energy, and we will get a stable OHe-nucleus bound state of length Rj, and 
energy ^(-Rf,). If there is no minimum, then we will have to conclude that 
no stable bound state can form at those distances. 

In the following, we shall go to 3rd-order perturbation theory. We as- 
sume that Hq has a spectrum \rp^ > of eigenfunctions with eigenvalues E®, 
and we assume that the unperturbed energy level is non-degenerate. The 
formulae for the wave function at order 2 and for the energy at order 3 are 
given by: 
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E n — E n 



+ <^\W\^> 



+ E 



<^ P;i \w\^ n > 



2 



E° - E° 

i,p^n n P 

< > £ ( go _U 2 ( 16 ) 

i,p^n V ,p'j^n 



IV'n > = |< > 



y <<«>, 

i,P7^n ^ n P' 



1 ^ l<<Wn>l 2 |)0 



i,p^ni' ,p'jtn (-^n ~~ -^p) ( -^n — -^V 



In our case, the non-degenerate unperturbed energy E® is the ground 
level of the hydrogenoid OHe atom: 

E° n = E OHe = -\m He (Z Z He a) 2 ~ -1.58 MeV (18) 

and the unperturbed eigenfunction \ipp i > are those of the hydrogen atom: 

i>l,M) = = RnM)Yr(0i,<pi), (19) 

where the Y™ are the normalised spherical harmonics, and where the radial 
part R n i is given by 

Rn,l{n) =C n ,l*r[ V Cg - e -o , (20) 



C n / being the normalization coefficient of R n i and ro being the Bohr radius 
of the OHe atom. The coefficients c q in (l20j) are recursively given by — = 

c/(< ? +2«+l) ' 
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R (fm) 

Figure 10: AEq = Eq(R) — EoHe (keV) up to orders 1 (upper), 2 (lower) 
and 3 (middle) for an external sodium nucleus, as a function of its distance 
R (fm). V = 30 MeV and a = 0.5 fm. 

4.1 Correction to the OHe energy 

First, we consider the effect of an approaching sodium nucleus on the OHe 
energy. Fig. [TU] shows the results for AE = E (R)-E O He for V = 30 MeV 
and a = 0.5 fm. We see that order 1 doesn't bring a large modification to 
the unperturbed energy (~ 10~ 3 keV for R between 50 and 15 fm), while 
order 2 gives the largest correction (~ 1 — 10 keV for R between 50 and 
15 fm). This change from order 1 to order 2 justifies the inclusion of order 
3 in the calculations, but it turns out that this one doesn't modify greatly 
the results from order 2, and that is why order 4 has not been added. We 
see on Fig. [10] that AEq is always decreasing, in other words that there 
is no minimum in this curve in the region of validity of the perturbative 
calculation. 

Similar results hold if one strengthens the nuclear potential, or if one 
considers different nuclei, as shown in Fig. [11] in the case of iodine. 

4.2 Interaction with the incoming nucleus and polarization 

In the same way, we can calculate the electrostatic and nuclear interaction 
energies between the perturbed charge distribution of the helium nucleus 
and the charge distribution of the A nucleus as a function of its distance R, 
as well as the mean position of He along the z— axis, that is the polarization 
of the OHe under the influence of the external nucleus. 

The electrostatic interaction energy between two charge distributions is 
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R (fm) 

Figure 11: AEq = Eq(R) — EoHe (keV) up to orders 1 (upper), 2 (lower) 
and 3 (middle) for an external iodine nucleus, as a function of its distance 
R (fm). V = 45 MeV and a = 0.5 fm. 



given by 

E el =[ [ K^^ drirt (21) 
JviJV2 \ri-r 2 \ 

where each integral is performed over the extension V\ or V2 of the corre- 
sponding charge distribution and where p\ and p 2 are the charge densities 
of each distribution. In the following, we take the first-order version of tp 
for He (as it is responsible for the dominant second-order shift in energy), 
and pi (ri) = \tpHe (jt) \ 2 , while the external nucleus is treated as a uniform 
sphere. 

The results are compared to the Stark potential used in the previous 
section in Fig. [12j We see, as might be expected, that the simplifying as- 
sumption of the constant electrical field for the nucleus is reasonable at large 
distance, while the gap becomes more pronounced around 15 fm, because 
the uniform Coulomb field is always stronger than the true one. The fact 
that E e i becomes repulsive at shorter distance is due to the change of the 
polarization of the OHe atom under the influence of the nuclear force of 
the sodium nucleus, which makes the helium component turning to positive 
mean Z\, that is, towards the external nucleus. 

We can also integrate the Woods-Saxon potential ,~ v ° r, over the 

fa ^ l+e( r ~ R A)/°- 

distribution of the helium nucleus to get the total nuclear interaction energy 
E nuc i- Adding it to the previous contribution gives us the curve of Fig. [T3l 
which has no sign of a potential barrier. 

Finally, we can calculate the mean value of the position z\ of the helium 
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Figure 12: Electrostatic interaction energy E e i (MeV) at order 1 (upper) 
compared to the Stark interaction energy Vstark (MeV) (lower) as a function 
of the distance R (fm) of an external sodium nucleus; Vo = 30 MeV and 
a = 0.5 fm. 
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Figure 13: Total interaction energy E e i + E nuc i (MeV) at order 1 as a 
function of the distance R (fm) of an external sodium nucleus; Vo = 30 MeV 
and a = 0.5 fm. 
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Figure 14: Polarization < z > (fm) at order 1 as a function of the distance 
R (fm) of an external sodium nucleus; Vb = 30 MeV and a = 0.5 fm. 

nucleus along the z— axis, that is the polarization of the OHe atom, which 
is simply obtained by 

< z > = < ll>He\zi\ll>He >= J M Z^lpHe (^l) ? (22) 

Fig. I14lrepresents the evolution of this polarization as a function of R. It can 
be seen that, at large distance, the polarization is negative due to Coulomb 
repulsion between nuclei, as expected. Thus, Fig. [14] shows that the OHe 
atom gets polarized, for R < 10 fm, in a direction that could allow repulsion, 
provided that the nuclear force is not already too strong at such distance. 
The addition of the nuclear interaction with Vb = 30 MeV and a = 0.5 fm 
in Fig. [131 shows that this condition is in fact not satisfied, giving rise to an 
attractive force at all distances. The modification of the nuclear parameters 
Vq and a (for example Vb = 10, 100, 200 MeV, a = 0.5, 1.5 fm), as well as 
the external nucleus, doesn't radically change the results, modifying only 
the distance from which the potential falls to nuclear values. 

5 Conclusion 

The advantages of the OHe composite-dark-matter scenario is that it is 
minimally related to the parameters of new physics and is dominantly based 
on the effects of known atomic and nuclear physics. However, the proper 
quantum treatment of this problem turns out to be rather complicated and 
involves several open questions. 
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We have presented here the state of the art of our studies of the nuclear 
physics of the OHe atoms, and found a difficulty in proving the original 
assumption of a potential barrier developing between OHe and the nucleus 
A, both classically and in perturbation theory. 

Open questions for further analysis nevertheless remain: 

(a) for the distances under consideration the size of the He nucleus may 
not be negligible and it may not be sufficient to treat it as a point-like 
particle; 

(b) beyond the nucleus the nuclear force falls down exponentially but it 
may be strong enough to cause a non-homogeneous perturbation of 
the OHe atomic ground state; 

(c) the nuclear force indeed leads to a change of the OHe polarization that 
might result in the creation of a dipole Coulomb barrier, as shown in 
the perturbative calculation, but this happens when the perturbative 
approach is no longer valid, and one should thus solve the Schrodinger 
equation numerically in this regime. 

The answer to these open questions may be crucial for asserting the nuclear- 
physics basis of the OHe model. If there is no dipole Coulomb barrier 
between OHe and nucleus, one gets a spectrum of states, which could have 
transitions to each other. Although the spectra we showed in the third 
section are not reliable in the nuclear region, it is clear that a particles will 
have nuclear bound states. Without a barrier, their transitions to them will 
be fast and dramatic. 

Hence, the model cannot work if no repulsive interaction appears at 
some distance between OHe and the nucleus, and the solution to the open 
questions of OHe nuclear physics is vital for the composite-dark-matter 
scenario. 
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